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We examine the dynamical property of the polarization function for electrons with the 
tilted Dirac cone, considering the case where the chemical potential measured from the 
contact point is finite, as found in the organic conductor a-(BEDT-TTF)2l3 (BEDT- 
TTF=bis(ethylene-dithio)tetrathiafulvalene) at ambient pressure. Using the tilted Weyl 
equation with the analytical treatment of the particle-hole excitation, we show that the 
polarization function as a function of both the frequency and the momentum exhibits cusps 
and nonmonotonic structures. The polarization function depends not only on the magnitude 
but also on the direction of the external momentum. These properties are characteristic of 
the tilted Dirac cone, and are in contrast to those in the isotropic case of graphene. Further- 
more, the results are applied to calculate the optical conductivity, plasma frequency, and the 
screening of Coulomb interaction, which are also strongly influenced by the tilted cone. 

KEYWORDS: Dirac cone, tilted Weyl equation, zerogap state, a-(BEDT-TTF)2l3, polarization 
function, optical conductivity, plasma frequency 



1. Introduction 

The recent discovery of graphene^^ has attracted much attention in the field of condensed 
matter because graphene exhibits massless Dirac fermions. Using the Weyl equation,^' 
which describes the motion of massless Dirac fermions, several anomalous electronic prop- 
erties have been investigated for a long time.^'^-* Another massless Dirac fermion is found in 
the quasi- two-dimensional organic conductor a-(BEDT-TTF)2l3 (BEDT-TTF=bis(ethylene- 
dithio)tetrathiafulvalene) under pressure.^' The existence was demonstrated theoretically®' 
using the band calculation, which is based on the transfer integrals estimated from the X-ray 
structure analysis. ^"^^ Such an energy band has been confirmed by first-principles calcula- 
tions.^^' This novel state elucidates a long standing problem of anomalous phenomena 
observed in the conductor under pressure. ^^'^^^ 

The massless Dirac fermion in the organic conductor is expected to exhibit a noticeable 
property due to the anisotropic velocity of the tilted Dirac cone, which is described by the 
tilted Weyl equation. ^^"^^^ It has been shown that the tilt affects the characteristic temper- 
ature dependence of the Hall coefficient^^^ where the magnitude of the tilt is characterized 
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by the parameter a = vq/vc {0 < ct < 1), the ratio of the tilting velocity to the cone veloc- 
ity. A recent experiment on the transport phenomena of a-(BEDT-TTF)2l3 suggested that 
a = 0.8 using theoretical evaluation. ^"^^ Since the conductor is a layered two-dimensional 
massless Dirac fermion system, ^^'^^^ the interplane magnetoresistance also exhibits noticeable 
properties, as shown theoretically by the angle dependence of the magnetic field. Moreover, 
new phenomena induced by the tilted Dirac cones have been maintained by calculating the 
transport coefficient under strong magnetic field, the electric-field-induced lifting of the valley 
degeneracy, ^"^^ and the easy-plane pseudo-spin ferromagnet leading to the Kosterlitz-Thouless 
transition. In addition to the above studies, dynamical properties such as electron-hole ex- 
citation and collective excitation are promizing ingredients for verifying the role of tilting. 
Although the electronic state has been studied extensively, the dynamical properties asso- 
ciated with the polarization function are not yet clear enough compared with those in the 
isotropic case of graphene. 

In the present paper, we analytically examine the dynamical polarization function with 
the arbitrary wavevector and frequency, and compare it with that in the isotropic case of 
graphene. ^^"^"^^ We examine the metallic state where the contact point of the Dirac cone is lo- 
cated below the Fermi energy as expected for the organic conductor a-(BEDT-TTF)2l3.^'^'^''^ 
In §2, formulation for the polarization function is given. In §3, the analytical expression of 
the imaginary part is calculated, while the real part is estimated semianalytically using the 
Kramers-Kronig relation. In terms of these results, we demonstrate the cusps and nonmono- 
tonic structures as characteristic of tilted Dirac cone. On the basis of the results of §3, we 
calculate the optical conductivity, plasma frequency, and screening of the Coulomb interaction 
in §4. The conclusion is given in §5. 

2. Formulation 



We consider the zero-gap state in a-(BEDT-TTF)2l3, which has two tilted Dirac cones. 
Among two contact points, ifcoi corresponding to two valleys of cones, we focus on one of 
them, which is given by the state located close to fco(= ^Ox) ^Oy) with kQx < and kQy > 0. 
For such a state, the effective Hamiltonian is expressed as^^-* 

77'=1,2 

which gives the tilted Weyl equation. The quantity a^^ denotes a creation operator of the 
electron where the momentum k is measured from that of the contact point. The apex of the 
Dirac cone for the conduction band touches that of the valence band. Using the Luttinger- 
Kohn representation, where the two states at the contact point are chosen as the basis with 
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Fig. 1. (a) Energy dispersion of the Dirac cone tilted along fca;-direction. The chemical potential is 
situated on the upper cone where the white circle denotes the Fermi surface. The angle Of, , which 
is measured from the k^-axis, is given by 9^ = aTcta,ii{ky / k^) in the k^-ky plane. (b)Chemical 
potential /i (mcV) as a function of the uniaxial pressure Pa (GPa) for a-(BEDT-TTF)2l3. 



the index 7, the matrix Hk is expressed as 

H,^k.( "» "'-'"^V (2) 

\Vi + IV2 f / 

Two states of the basis are chosen so as to give vq = (vq,0),vi = {v,0),V2 = {0,v') with 
V = v' = Vc (we use h = 1). DiagonaHzing eq. (2), the energy around the contact point is 
obtained as 

= vok^ + svc\k\, (s = ±l), (3) 

where the first term comes from the tilting along the /c^^-direction. Defining a tilting parameter 
a{< 1) as 

eq. (3) is rewritten as 

Lk = {s + a cos 6k)vck, (5) 

where k = {kx,ky), k = \k\ and tan^^ = ky/kx- From eq. (5), the trajectory for the fixed 
^s,k gives the ellipse with a focus at A; = 0. Figure 1(a) shows the energy dispersion obtained 
from eq. (3) where the Dirac cone of the upper band touches that of the lower band at the 
contact point, = 0. In the present paper, the energy is scaled by the Fermi energy (i.e., the 
chemical potential /i), which is measured from the contact point. In Fig. 1(b), ^"^^ the chemical 
potential of a-(BEDT-TTF)2l3 is shown as a function of the uniaxial pressure Pa along the 
a-axis. With increasing Pa from ambient pressure, fx decreases from 7 meV and reduces to 
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zero at Pa = 2.9 kbar, above which the zero gap state is obtained. The case of ^ = is briefly 
mentioned in the discussion. The eigenvector for eq. (3) is given by 

^,W^^(J. (6) 

Thus, eq. (1) is rewritten as 

X] (^fc)77' = X] is,kw\ ,^Ws^k, (7) 

77' = 1,2 s=± 

where 

= ^(F,(fc))^a^,fe. (8) 

k 

For these two operators, Green's functions of the single particle are defined by 

»1/T 



GsiKiEn) = -\ I dT (TrWs,k{r)wsMO)^) e"""", (9a) 

G};^,{k,ien) = --J AT(Tra^^k{T)a^',kW) (9b) 

where en{= i(2n + 1)ttT) is the Matsubara frequency for the fermion. T is the temperature 
and Tr is the ordering operator for the imaginary time (we use /cb = 1) • Equation (9b) is 
written explicitly as 

The polarization function per valley is calculated as 

-l/T 



U{q,iu^) = - 2^/ '^^e— (T.p^,,(r)p^, ,^(0)> 
2 

2 



77 

2 

— Yl ^^n' (^ + 9' + iwm)G^,^(/2, ien) 

kwyy' 

— ^ |F;(fc)F,,(A; + q)|^G,/(fc + q,ie„ + itJ„)G,(fc,ie„), (11) 

knss' 

where the freedom of the spin is included. ia;m(= \2rmiT) is the Matsubara frequency for the 
boson, p^^q is the electron density operator defined as p^^q = "Y^j^ a)^ j^a^^k+q and 

\F^ik)F,,{k')f = \[l + ss' cos(0fc - 9k')] . (12) 
After performing an analytical continuation given by icom. — > w + ir], (r] = +0), one obtains 

\tpUi,\TP (I. , f{is,k) - f{is',k+q) 



1^^-^ ' k ' UJ + IT)- {is',k+q - is,k) 

= Y,^sAq,^). (13) 
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where /(^) = (1 + e^^~^^^'^^)~^ . The term Ilgs comes from the electron-hole excitation of the 
intraband, while the term H^^-s comes from that of the interband. 

In the present paper, we examine the case for T = and calculate eq. (13) analytically by 
performing the integration over k. The case for u > and /_f > is calculated, while the case 
for < is obtained from 

U{q,-u)=[U{-q,uj)]*, (14) 

which comes from the property of py^q of the tilted cone, i.e., without inversion symmetry. In 
the isotropic case with inversion symmetry, one obtains n(q, —uj) = [n(q,a;)]*. 

Here, we note that, by taking account of the freedom of both spin and valley, the polar- 
ization function of the total system is given by 

U'-'-\q,oj)=Uiq,oj)+Ui-q,uj). (15) 

The quantity of I[{q,uj) per spin and valley is calculated in §3, while charge response is 
calculated using n*°*^'(q, tj) in §4. 

3. Polarization Function 

3.1 Analytical calculation 

We examine the imaginary part using eq. (13). By noting Imll = Imll^ = due to the 

valence band fully occupied, we calculate the remaining parts of Imn_|„+ and Imll |_, which 

correspond to the process excited from the conduction band and valence band, respectively. 
After performing the tedious but straightforward calculation, we obtain the analytical results 
(Appendices A and B): 

lmU++{q,eq,u;)= (16a) 

C=1A,2A,3A 

lmU^+{q,eq,uj) = U'l. (16b) 

C=1B,2B,3B 

The results consist of six regions on the q-oj plane where the typical case with 9g = 7r/2 and 
a = 0.8 is shown in Fig. 2. These regions of lA, 2A, 3A, IB, 2B, and 3B are classified into 
two regions, A and B, corresponding to the process of intraband and interband excitations, 
respectively. The regions A and B are separated by a solid line expressed as 

Wres = (1 + a cos Oq)Vcq, (17) 

which is called the resonance frequency. The resonance frequency is obtained owing to the 
nesting of the excitations with the linear dispersion, and the polarization function diverges 
with the chirality factor eq.(12) taking a maximum. '^^^ The boundary between 2A and 3A is 
given by uj+. The boundary between 2A and lA ( 2B and IB) is given by coa- The boundary 
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Fig. 2. Several regions on the q-u) plane for Imn(q,a;) with the fixed 9q = 7r/2 and a = 0.8. The 
solid line denotes the resonance frequency, lo = (^vcs- The dash-dotted line between 2B and 3B 
corresponds to wb and the line between lA and 2A (and also fB and 2B) corresponds to oja- The 
imaginary part is absent in the regions of both 3A and IB. The boundary between 2A and 3A is 
given by cj+ . The intersection point of uja and u-^cs is given hy Vcq/ ^J■ ~ 1/(1 + Q! cos 9q ) and uj/ fi ~ 1 
while wa/a^ — 2/(1 + a) at Vcq/ H = 0, and = at v^qj ~ 2/ (1 — a^)(l — o? cos^ dq). All 
the boundaries except for the solid line are followed by the cups, i.e, the jump of the derivative 
with respect to q and w. 



between 2B and 3B is given by wb- These frequencies are calculated as 



±- = a COS^q - 2 + V r 2 + 1 2 ' (1^) 



In the case of the isotropic Dirac cone (i.e., a = 0),28'29) there is a boundary given by 
u}/fi + Vcq/fJ- = 2 which separates lA and 2A (IB and 2B) for the intraband (the interband). 
In the regions 3A and IB, the imaginary part vanishes. In the present case of the tilted Dirac 
cone, the boundary between lA and 2A exhibits a noticeable behavior characterized by the 
appearance of cusps for the imaginary part as shown in this section. 

The real part is calculated from the imaginary part using the Kramers-Kronig relation: 

1 /-^ lmUiq,9^,x) 



1 f°° 
ReU{q,eq,U}) =-V dx 

J-oo 
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Fig. 3. (Color online) Normalized imaginary part, Imll{dq, uj)v^ / fi, as a function of uj/fi with 
q = 2/i/uc and 6q = it/2. This behavior is expected in the case of 1/(1 + acosdq) < Vcq/h < 
2/ {y/l — y/l — cos^ 6q). The region of < w < Wios (i^ > i^rcs) comes from the intraband 
(interband) excitation. The dotted line represents the isotropic case. The cusp in the intraband 
region is found at uj = uja, while that in the interband region is found sd, lu = ujb (the inset). 



■■""fa."..'-) ^ WJ, ■mn(,., + (>..x) 

TTJo X - (jj VrJo X + UJ 

where we used, from eq. (14), the relation 

Ren(g, —uj) = Rell{—q,uj) = Ren((/,7r + 9q,uj), 

Imn(q, -Lo) = -ImU{-q,uj) = -Imn(Q',7r + 6g,uj). (21) 
The result of the semianalytical calculation is given in Appendix C 

3.2 Behavior of the polarization function 

The imaginary part exhibits a variety of uj dependences depending on the magnitude of 
q, which is divided into the following three regions. In the region of the small q, a gap exists 
in the interband excitation, while in the region of the large q a gap exists in the intraband 
excitation. The intermediate region of q is located between these two regions. 

First, we show a novel behavior seen for the intermediate magnitude of q, which comes 
from the tilted Dirac cone. In Fig. 3, lmll{6q,uj)v'^/ with fixed v^q/^i = 2 and 6q = 7r/2 
is shown as a function of uj/fi. There are two cusps at a; = ua in the intraband region 
(0 < a; < uJrcs) and at a; = in the interband region (w > Wrcs), where uja and ujb are given 
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by eq. (B-13). Note that 

< ^rcs < ^B, (22) 

since regions for intraband and interband excitations are separated by Wres- 

Here, using Fig. 4, we explain the cusp at w = oja in the intraband region. Equation (13) 
shows the following three kinds of conditions for obtaining the imaginary part. The first one 
is the energy conservation given by 

^ = i+,k+q - i+,k, (23) 

and the second one is the possible process of the electron-hole excitation given by f{^+^k) ~ 
f{(,+^k+q) = li which leads to 

e+,fc < (24) 

(+,k+q > ^J'■ (25) 

The third one, which is discussed later is a factor given by eq. (12). From the second condition, 
it is found that k is allowed inside of ellipse (I) and outside of ellipse (II) on the kx-ky plane 
where ellipse (I) and ellipse (II) denote E,+,k = /U and ^+,fc+q = ^, respectively. By defining 
R as the intersection point between ellipse (I) and ky = kx tan 9q{> 0), one finds oj = uj^es 
being the maximum on the line given by ky = kxtan9q{> 0) (i.e., on the line OR). The cusp 
is understood as follows when the origin O (i.e., kx = and ky = ) is located outside ellipse 
(II). We define Pi and P2 as the intersection points between ellipse (I) and ellipse (II) at which 
u{= k — C+,k+q) 0. Thus, uj increases from zero to cOrcs when the point moves from Pi{ 
or P2) to R on ellipse (I). Here, we consider ellipse (III), which has the same focus as ellipse 
(I) and touches ellipse (II) at the point A. It can be shown that such ellipse (III) is given by 

= fi-UA, (26) 

where 0JAi> 0) is equal to that found in Fig. 3. By noting ^+,fc+<3 — /i > at fe = 0, and the 
existence of Pi and P2, the condition for the cusp at ua is given by 

^ < v,q/n < , (27) 



l + acos^q "'^ ^(1 - a2)(l — cos^ 6„ I 
where Wres > > 0. On ellipse (II), u takes a minimum at Pi (and P2) and a maximum 
at A suggesting a saddle point at A, which gives rise to the cusp. Actually, this leads to an 
excess contribution for the imaginary part, which becomes singular for uj — >• loa — 0, i.e., 

6 (lmn{q,9g,uj)) oc -y/u)A - (28) 

Next we examine the imaginary part in the other cases of the small q and the large 
q, which are found for Vcq/fi < 1/(1 + acos6'q) and 2/(Vl — a^\/l — cos^ 9q) < Vcq/fi, 
respectively. These examples are shown in Figs. 5(a) and 5(b) by choosing q = fj,/2vc (a) 
(upper panel) and q = 4^1 /vc (b) (lower panel). Note that the cusp is also found at uja or 
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Fig. 4. (Color online) Fermi surfaces for ^fc = ^ (I), Cfc+g = M (H), and = fJ' — i^A (M) on the kx-ky 
plane. Point A where I and II contact gives rise to the emergence of the cusp. Pi and P2 are the 
intersection points between ellipse (I) and ellipse (II), and O is the origin. The condition for the 
resonance is obtained for ky = kx tan^g, which is found on the line connecting O and R. 



a;+, which can be understood from the saddle point similar to that in Fig. 4. There is a 
broad peak for u < cjrcs, which is also characteristic of a tilted cone and is absent in the 
isotropic case. The location of such a peak is almost proportional to q. A noticeable difference 
is seen between the small q and the large q (i.e., between the case q = fj,/2vc and that of 
q = Afji/vc)- For q = fi/2vc (upper panel), the interband contribution of lmn(0q, a;)t;^//i is 
very small and increases from zero with increasing uj{> oja)- As seen from Figs. 3 and 5 (a), 
Imll(9q,uj) of the intraband excitation (a; < oj^cs) is much larger than that of the interband 
excitation (cj > Wrcs)- An opposite behavior emerges when Vcq/fJ- >> 1 (Fig. 5 (b)). For 
q = 4:fi/vc (lower panel), the interband contribution of Imn(^q, uj)v'^/fi is large and is close to 
u ~ Wrcs- Compared with that in the isotropic case, the intraband contribution is enhanced 
but the interband contribution is suppressed. Such a fact may be ascribed to a sum rule 
given by dujujlmll{q,u}) = (l/8)g'^A, where A is the upper cutoff of the frequency. For 
Vcq/^-i < 1/(1 + acos^q), where uja > f^resi the imaginary part is absent in the interval range 
of 

Wjcs < < i-^A- ^otl a small cu, the imaginary part is convex downward in the isotropic case 
owing to the absence of the backward scattering, while it is convex upward in the tilted case 
owing to Pi and P2 located away from the backward scattering. Note that the cusp at a;^ is 
present for the arbitrary q, as seen from Fig. 3. 
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Fig. 5. (Color online) Normalized imaginary part, lmll{0q,uj)v^/ ^, as a function of w/zi in the case 
oi 9q = tt/2 for q = i-i/2vc (a) (upper panel) and for q = (b) (lower panel). The behaviors of 

panels (a) and (b) are found for Vcq/n < 1/(1+Q:cos0q) and — cp-yj\ — cP- cos^ 0q) < v^qj [i, 

respectively. The result in the isotropic case is shown by the dashed curve. 



We examine the 0q dependence of the imaginary part. On the basis of eq. (27), the region 
for the existence of uja is shown on the v^qj^i-Oq plane in Fig. 6. The behavior for the fixed 
^q(= ^/2) is as follows. In the region for the small Vc(?//i, the cusp at a; = liJa appears in 
the interband process (Fig. 5(a)), while in the middle region with an intermediate Vcqj^i^ oja 
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exists in the intraband process (Fig. 3). For the large Uc^Z/i, the cusp appears at u^, which 
is the lower bound of the imaginary part (Fig. 5 (b)). Note that the cusp is found, except for 
9q = and vr. 

Here, we briefly mention the 0q dependences on oja, ^tcs, and u)b with a fixed Vcq/fi- 
With increasing 6q from zero to vr, both lob and Uj-^s decrease monotonically, but up^ remains 
almost constant until uja = i^rcs- The cusp at ua moves from the intraband region to the 
interband region, where the boundary is given by oj^cs- Those frequencies also depend on the 
degree of tilting, a. With increasing a from zero to 1, both and increase monotonically. 
The resonance frequency Urcs increases (decreases) for < 0q < 7r/2 (7r/2 < < vr) with 
increasing a, while oj^cs at ^q = vr/2 is independent of a. 

Figure 7 shows the normalized imaginary part lmll{6q, u})v'^ / fi on the plane of Vcq/fJ- and 
u/n for 6q = 0(a), vr/2(b), and vr (c). The color gauge with the gradation represents the 
magnitude of the imaginary part. The global feature is mainly determined by the property 
of 

Wres- In the case of VcQ/ ^ >> 1 where the characteristic energy becomes much larger than 
the interband energy (~ /z), lmn(0q,(x>) of the intraband excitation {oo < oJj-cs) becomes much 
smaller than that of the interband excitation (u > Wrcs) in contrast to the case where Vcq/fJ^ = 2 
(Fig. 3). The broad peak in the intraband excitation (i.e., uj < Wrcs) does not change much for 
9q < vr/2, while it is strongly masked for vr/2 < 0^ < vr due to the rapid decreasing lOj-qs. The 
cusp at w = UA is seen for < < vr, as seen also from the analytical condition of eq. (27). 
The cusp in the interband process does exist in an intermediate region of q for the arbitrary 
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9q (except for and vr). 

Finally in this subsection, we examine the real part, which is calculated using eq. (20) 
(Appendix C). The numerical results of Kell{q,6q,uj)v'^/ fj, corresponding to the imaginary 
part of Figs. 3, 5(a), and 5(b) are shown in Fig. 8. When q > {fi/ Vc) / {1 + a cos 6q) as seen from 
Figs. 8(a) and 8(c), Rell{q,9q,uj) diverges for cj — curcs — 0. The cusp found in the imaginary 
part also appears in the real part, e.g., wa and ojb for Vcq/h = 2, ujb for Vcq/fJ- = 1/2, and 
a;+ and lo-q for Vcq/fJ- = 4. A noticeable difference compared with the isotropic case is the 
nonmonotonic behavior for cj < Wrcs and Vcq/fJ' > 2/(1 + a cos 9q), e.g., exhibiting a minimum 
for Vcq/fJ^ = 2 and 4. In the regime of the interband { uj > lo^cs), the real part is similar but is 
slightly small compared with that in the isotropic case. 

4. Optical Conductivity, Plasma Mode, and Screening Effect 

4-1 Optical conductivity 

The optical conductivity is obtained from two kinds of electron-hole excitations, i.e., the 
intraband processes and interband processes. The intraband excitation gives a Drude-like 
behavior in the presence of the impurity scattering, while the interband excitation gives a 
gaplike behavior. In this subsection, we calculate the interband optical conductivity, while the 
effect of the intraband one is left for discussion. 

In terms of the imaginary part of n(q,w), the optical conductivity cr{uj) is calculated as^^^ 

• 2 

tC CO 

a{io) = lim—^{U{q,u)+U{-q,u;)}, (29) 

where the RPA gives the same result. Since lmn(0,(x>) becomes finite only for the interband 
electron-hole excitation, Jlea{uj) is calculated as (Appendix ) 

Reo-(w) =0- G l^^-u 



1 + a 



e2 



87rh 



" + ' l + a|cosg,| - "J ^<^'^ ) + ''''' [ l-a\coseq\ " " ' ^<^-'^^ 



1 — Q / \ 1 + a 
where G is a step function, G<(x) = xy/ 1 — x"^ — arccos(3;) for |x| < 1, and are defined 

by 




(31) 



In Fig. 9, the normalized optical conductivity Kea{6q,Lo)Ah/e'^ is shown as a function of Lo/fi 
for 6q = 0, 7r/4 and tt/2, respectively, with a = 0.8. The dotted line corresponds to the 
isotropic case where the normalized Ilea{9q,uj) is zero for u < fi and is constant for fi < oj. 
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The conductivity I{ea{6q,u}) is zero for lo < 2^/(1 + q) and 1 for uj > 2/i/(l — a), while it 
takes an intermediate value for 2^/(1 + a) < w < 2^/(1 — a). The reason for such a behavior 
is illustrated in Figs. 10(a) and 10(b). The electron-hole excitation for the conductivity is a 
vertical transition from the valence band to the conduction band due to g = 0. Figure 10(a) 
denotes the process with a fixed energy oj = — The contour projected on the k^-ky 
plane is shown by a circle in Fig. 10(b). The ellipse of the thick line denotes the Fermi surface 
with ^ = ^. The electron-hole excitation is allowed outside of the thick line since the hole 
is created above the Fermi surface of the conduction band. Thus, the process is completely 
forbidden for uj < 2fi/{l + q), while such a process occurs on all the corresponding circles for 
Lo > 2/i/(l — a). In the intermediate case of 2^/(1 + a) < u < 2/i/(l — a), the electron-hole 
excitation occurs on the partial part of the circles outside the thick line, i.e., for 9q satisfying 
the condition u > 2/i/(l -|- acos^q). The tilting of the Dirac cone is essential for the optical 
conductivity that depends on 6q . 

The conductivity can be calculated directly by taking the limit of g — )■ in the imaginary 
part of the dielectric function. From eq. (13), the polarization function, which includes the 
excitation from the valence band to the conduction band is written as 

(27r)2y \ \k\\k + q\j u + ir] + ^^^k^q-^+^k ' 
By expanding up to 0{q'^), the imaginary part for the small |q| is rewritten as 

[a I'aui 1 -cos2(93-0q) 2 



Imn_+(g,^) ^ jd^jdkk q 

X 5 [ 2k — aqx — q cos((p — 6„) \ Q { k + ak cos ip — — 



vr 



a; {2tt)'' 



^2 



J ^ \ 1 + a COS if ) 



=|^i^c(^J. (33) 
The quantity Fc{q) is zero for u < 2fi/{l + a) and Fc{q) = 1 for 2;u/(l — q) < cj, and 

F,{dq) =- {1 - cos\y^ - 6q)} 

= ~ (^^0 ~ \ '^i^ 2v3o cos 29^ (34) 



for 2/x/ {1 + a) < oj < 2/z/ (1 — a), where 

T = 



cos¥Po = -f— -1). (35) 



Thus, we obtain the optical conductivity as 



e2 



^(^) = ^Pc{eq). (36) 

One finds the following 9q dependence of the optical conductivity in the interval region of 
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2;u/(l + a) < a; < 2/i/(l — a). From eq. (34), a with a fixed uj takes a maximum (minimum) 
at 0q = 7r/2 and a minimum (maximum) at 9q = and n for u < 2fi {co > 2fi). 



4-2 Plasma mode 



We examine the plasma frequency by adding the Coulomb interaction. Within the RPA, 
the plasma frequency can be calculated from the pole of the polarization function, which is 
obtained by 



l + VqReU{q,eq,u}) = 0, 

where Vq is the Fourier transform of the Coulomb interaction, 

27re2 



(37) 



(38) 



Equation (37) is written explicitly as 
2 



1 



-V, 



1^2 Yl 

k,s=+,s'- 



F^{k)FAk + q] 



I 2 /(6fc) - /(6'fc+q) 

UJ + il]- (6'fc+g - i 



sk 



1 



-V, 



"^(2^)2 

k{k 







2lT 



dkk dipfi^u) 



1 + 



1 



1 + 



k{k 



\k\\k 



1 



ik + ik+q 



(39) 



\k\\k + q\] u} + ik- ik+q 
which is calculated using the real part of n(g, ^q,a;) obtained in §3. In Fig. 11, plasma fre- 
quency as a function of v^qj is shown with several choices of Oq. In order to understand the 
behavior with a small we calculate the collective mode by expanding eq. (39) in terms of 
q. The lowest order is obtained as 

2li^ 



eo 



-C{9q)\q\ 



(40) 



where the coefficient C{9q) is given by 



C{9, 



1 



2-7T 

- I dip 

vr 



1 + ^ 



1 - COS^((^ - 9q) 

1 + a cos p 
1 + 



a sm 



Vl -0^(1 + VT 



(41) 



a^[i -\- V L — ^ 

In this case, the plasma frequency takes a minimum, 2/(1 + \/l — a^), at 9q = and vr, and a 
maximum , 2/(\/l — 0^(1 + \/l — o^)), at 9q = tt/2. The plasma frequency in Fig. 11 reduces 
to the result in eq. (41) in the limit of a small q. The global behavior including a large q is 
reported in a separate paper. 
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4-3 Screening of Coulomb interaction 

We examine the static screening in the presence of Coulomb interaction. In the static case, 
i.e., uj = 0, the effective interaction within the RPA is given by 

v^^\q,e„0)= 7 (42) 

1 + Vqll{q,eq,0) 

The real part of Il{q, Og, 0) is shown in Fig. 12 with the fixed 9q = 0, 7r/4 and tt/2 where the 
isotropic case is also shown by the dotted line. There is a cusp at 2kF(9q) (except for 6q = 0), 
which is given by 

kpiOn) = (43) 
^ ^(l-a2)(i_a2 cos2 9q) ' 

For q < 2k-p{0q), n(q,^q,0) remains constant as found in the conventional two-dimensional 
electron gas. In this regime, the interaction is screened as 

^;^P^(g, Oq, 0) ^ 27reV{eo('Z + ^tf)}, (44) 

where qtf is the Thomas-Fermi screening constant given by 2eV{eo(l - a^f/'^]. For a large 
q, we obtain 

n(g,g„0)^— = \ , (45) 

81^1 — COS"' Oq 

which corresponds to ;U = in eq. (C-2a). In the case of a large q, the property of the Coulomb 
interaction remains with the the dielectric constant replaced by 

7rp2 

eo + , (46) 

4t;c Y^l — q2 cos2 9q 

The dielectric constant takes a maximum (minimum ) at 9q = 0(7r/2) suggesting a fact that 
the anisotropy of the velocity gives rise to the enhancement of the charge response. Thus, we 
found that the effect of the tilted cone also emerges in the 9q dependence of the dielectric 
constant. 

5. Conclusions 

We have examined the property of the polarization function of the massless Dirac parti- 
cle with the tilted cone and a finite doping. Using the tilted Weyl equation, the dynamical 
polarization function with the external momentum q and frequency a; is treated analytically, 
and is applied to calculate the optical conductivity, plasma frequency, and screening of the 
Coulomb interaction. The tilting of the Dirac cone gives the following effects on the polariza- 
tion function, which is determined by the intraband and interband excitations. 

Several cusps as a function of the frequency, w, exist in both the imaginary part and 
real part of the polarization function at the characteristic frequencies, w = w^, ^b, and a;+, 
depending on the momentum q and the tilting parameter a. Nonmonotonic structures and 
cusps originate from the tilting and are in contrast to those in the isotropic case (a = 0). The 
cusps are understood in terms of the saddle point of the particle-hole excitation energy on 
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the plane of the momentum. The tilting effect also appears in the anisotropy of the resonance 
frequency cOres, which separates the intraband excitation from that of interband one. The 
intensity of Imn(^q,a;) exhibits a strong asymmetry between the intraband region and the 
interband region. In the case of Vcq/^ < 1, linll{9q,uj) of the intraband excitation (w < Wres) 
is much larger than that of the interband excitation (oj > lOj-qs), while the behavior is opposite 
for Vcq/fJ^ » 1. 

The optical conductivity, which is finite above the critical frequency increases continu- 
ously from zero and reaches a universal value. The plasma frequency for a small q, which is 
proportional to q^^'^, depends on 9q, and is symmetric with respect to 9q = 7r/2, where the 
coefficient takes a maximum at 9q = -k/2. Such an anisotropy also exists in the screening of the 
Coulomb interaction, where the renormalized dielectric constant becomes largest at 9q = 0. 
The effect of interaction on two tilted cones is reported in a separate paper within RPA. 

We comment on the characteristic energy relevant to a-(BEDT-TTF)2l3. It is expected 
that several cusps obtained in the present paper occur for the energy around the Fermi energy 
being ~ 5 meV. However, the observation of these cusps is not straightforward since the 
neutron scattering is inapplicable owing to the absence of magnetic ordering and the energy 
scale for the X-ray experiment is too high. One possible way is to use the electromagnetic wave 
between microwaves and infrared waves. The optical conductivity (Fig. 9), which is obtained 
in the zero limit of momentum, also exhibits cusps which are different from those of w^, ojb, 
and For the frequency dependence of the optical conductivity, the Drude like behavior of 
the intraband optical conductivity is well separated from the interband one since the energy 
for the impurity scattering (~ 0.1 meV) is much smaller than the Fermi energy (~ 5 meV). 
This is in contrast to the case of the graphene where the Drude like contribution becomes large 
in addition to the interband one^^^ owing to the large scattering energy of 10 meV.35) Thus, 
we may find the behavior showing that the optical conductivity increases continuously from 
the lower edge of hw = 2^/(1 + a) ( ~ 5 meV) and that the steep increase of the conductivity 
at the edge depends on 9q{^ 0). In addition to such a continuous variation, we may expect 
the effect of finite temperature. Although it is similar to Fig. 9, it can be distinguished from it 
since it approaches tangentially to the limiting value. Furthermore, there are also contributions 
from other two bands located around ~ 500 meV below the Fermi energy where a similar 
variation with respect to the frequency is expected. 

Finally, we briefiy mention the case of the zero gap state under pressure in which the 
Fermi energy is located on the contact point of the cones (i.e., /i — )■ 0). Although the effect of 
tilting still existsts, through the anisotropy of uj^cs = (1 + acos^q), Il{q,u}) in Fig. 2 remains 
finite only in a single region of the interband, i.e., lob moves to uj^es- The optical conductivity 
becomes constant for both the arbitrary u and 9q, while the plasma mode vanishes owing to 
the absence of the Fermi surface, as seen from eq. (40). 



16/28 



J. Phys. Soc. Jpn. 
Acknowledgment 



Full Paper 



The authors are thankful to R. Roldan, J.-N. Fuchs, M. O. Goerbig, F. Piechon, and 
G. Montambaux for fruitful discussions. Y.S. is indebted to the Daiko foundation for the 
financial aid to the present work. This work was financially supported in part by a Grant-in- 
Aid for Special Coordination Funds for Promoting Science and Technology (SCF), Scientific 
Research on Innovative Areas 20110002, and Scientific Research 19740205 from the Ministry 
of Education, Culture, Sports, Science, and Technology in Japan. 

Appendix A: Calculation of the Imaginary Part 

By noting that f{^s,k) ©(^ — for T — )• and using (x + = V\ — iTTS{x), the 

imaginary part of eq. (13) is rewritten as 



ImU{q,uj) = ^n^^,' 

ss' 

--^ jdk\Fl{k)F+{k + q)|' [g(/x - C+,k) - e(/i - 



-i- ldk\Fl{k)F^{k + q)f e(/x-^_,fcj 



+ ^ I dk\Flik)F^{k + q)f Q{^i-C+,k) 



6{UJ + - C+,k+q) 

e(/i - C~,k+q)\ S{^ + C-,k- ^^,k+q) 
0(/^ - ^-,k+q) S{UJ + - ^~,fc+q) 



+ i_ ldk\Fi{k)F+{k + q)f Gai-^-,fcj 



0(/^ - ^+,k+q)\S{uJ + - ^+,k+q), (A-1) 

where Imll = and Imll^ = owing to /x > 0. + q —7-= k' in Inin_|-+ and Imll |-. By 

making use of eqs. (3) and (12), eq. (A-1) is written as 

imn. 



2ttVcJ 2 

fdk'- 

27rvr.. 2L 



imn_ 



-.^[dk'J 

2-KVr.. 2 



1 + COs(6'fc - Ok+q) 
1 + COs(6'fe - 6'fc_q) 

1 - COs(6'fc - Ou-q) 



Vr. \ Vr, 



e 



k — akr 1 5 



Vr. 



Using the relation 



cos(6'fc - ek±q) 



VcJ V 



k ±q cos(6'fc — Og 



k+\k- q\ , 



k-\k 



(A-2a) 



k-{k±q) _ 
k\k±q\ ~ + g2 ± 2kq cos t/i. — d„] 



(A-2b) 
(A-3) 



with 0k (Oq) being an angle between k (q) and the x-axis, eq. (A-2) is rewritten as 



Imn 



++ 



27rt;, 



oo pi 

kdk 



cJo 



du 



X e 



ix/r= 

k — ak 



1 + 



k + qu 



y/W+(f + 2kqu 



ucosl 



v? sin ( 
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2'KV, 



kdk [ 



cJO 



1 + 



k — qu 



-y/fc^ + g2 — 2kqu 



X 9| k — ak 

Vr 



ncosf 



sm ( 



(A4a) 



imn_ 



2ttv, 



kdk i 



k — qu 



\/k'^ -\- q"^ — 2kqu 



X Q[ k + ak 



ncos 



9^- VT 



sin ( 



where (f = Ok — (^q, ^ = — = ^ ~ ck'Z cos Bq^ where gi,g2, and are given by 



gi{u) = u + k — \/k^ + q"^ + 2kqu , 



g2{u) = V — k + a/Zc^ + g2 — 2kqu , 



(A4b) 

(A-5a) 
(A-5b) 
(A-5c) 



gsiu) = V — k — _|_ g2 _ 2kqu . 
By noting the necessary condition for gi{u) = 0(— l<n<l), the 5 function is evaluated as 



gi{u)=Q^Q{q-u)Q\k 
g^iu) = ^ Qiq - iy)Q ( k 



q — u 
2 

q + v 



g^^u) = ^ e(z. - g)G [k - ^) e - A;) . 

The integration with respect to u is performed as 



(A-6) 
(A-7) 
(A-8) 



imn. 



2ttVc 



dk\ 



{2k + vf-q^ 



q-u 



^ 2fi + uj . ^ .2k + u 
X fc) [q + av cos Oq) 



+ Q sin 9q sj q^ — v'^\ 



2k 

Q 



2'KV, 



dk\ 



{2k -v) 



cJO 



2fi — uj 



q'^ — \ 2 



2k -u 



+ v 



X B I — {q + aucosOq) h Qsin^^y — i/ 



2 _ ,.2, / / ^ ^ 1 _ I 



(A-9a) 



imn_ 



1 rg!f-(2k-^) 



2'KV, 



,2 _ r,2 



W-,)e(.-ii^)e(^ 



X G 



2[i — ijj 



+ (g + az^ cos Bq 



2k -V 



asm( 



2k -V 



(A-9b) 
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imn_ 



47ruc y'q^ - 


9 




OJ 






1 q^Q{q 


-u) 




-1/2 




UJ 






1 q'^Q{iy 


-q) 



JdkWk'^ - le (A:' - 1) 

(g + az/ cos 6q) k' + a sin 0g \/ q"^ — k"^ — 1^ 

Jdk'vW^e {k' - 1) 

(g + az^ cos 0q) k' + a sin 0q \/ q"^ — k"^ — 1^ 



(A- 10a) 



c yJV^ - q^ 



dk'Vl - k'^e {k' + 1) e (1 - A;') 



X G 



2fl — Ld 



+ {q + avcos 6q) k' — asm.9q\/ — q^\J\ — k"^ 



(A-lOb) 



The step function in the above equation for ((7 > 0, > 0) gives the fohowing condition. 

When q> u, q+av cos > for < 0q < vr. When v > q, q+av cos > for < 0q < | 
, and q + aucosOq < with 00 > — ^"^^JlasT ^''^ § < 0q < vr. Thus, eqs. (A-lOa) and 
(A- 10b) are rewritten as 

1 q'^e{q-i^) 



Im n_i 



+ 



imn_ 



+ 



+ 



+ 



+ 





zy2 


1 q^e{q- 






Zy2 


1 q^e{q - 




4vrwc v^g2_ 


Zy2 


1 g20(^_ 




47r?;c y/iy2_ 




1 g20(j, _ 


-q) 






1 q'^e{u - 


-q) 


4:7rvc ^v'^- 


q' 


1 (^Q{v - 


-q) 


4:7rvc Y^zy2_ 


q2 


1 q^e{u - 


-q) 



■e{xi - i)Jdk'Vk'^ - 1 



e{x2 - l)jdkWk'^ - 1 



e 



OJ 



Vu)e{l- X2) / dk'Vk'^ - 1 



(A- 11a) 



6(1 - \x2\)@iq + aucosOq) / dk'^l - k"^ 

Jx+ 

Q{1 - \x2\)Q{-q - av cosOq) jdk'^/l-k"^ 



e [Vu 
e I Vu ■ 



— ^) e(-l - 2:2)9(9 + 01/ COS 0g) fdk'Vl-k'^ 
^ e(x2- l)e(-g-az/cos0g) jjk'Vl - k'^ 
Q{-l-X2)Q{q + avcoseq) jdk'^/l - k'^ 
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eiVu - - — ^ I e(x2 - l)G(-g - au cos 6^) I dk'Vl- k'^ 



, 1 q^e{iy-q) ^ fu-2fi 



(A-llb) 



Appendix B: Expression of the Imaginary Part 

Performing the ^'-integration in eqs. (A-lla) and (A-llb), eqs. (16a) and (16b) are calcu- 
lated as follows: 



n'lA = f{q,^) G>{xf) -hsgn(a;i - l)Gy{x^ ) - G>(x^) - sgn(x2 - l)G'>(a::2 ) Q{uJrcs - uj)@{uja 

(B-la) 

n'aA = /((?,i^)[g>(x+) + sgn(xi - l)G>(xr)]e(L^rcs -^^)e(L^ -L^A)e(a; -0J+), (B-lb) 
U'^j, = 0-eiu+-io), (B-lc) 



n'/g = O-e(w-a;rcs)0(^^A-w), (B-2a) 

^2B = fil^'^) + sgn(x2 - 1)G<(2;^) sgn(x2 + 1)G<(-X2 ) 0(a; - tJrcs)0(^^B - io)Q{uj - to a) 

(B-2b) 

n3B = /('7,Z^)2^0(^-^B). 



hi the above equations, we define 

1 



8ttVc ^|g2 _z,2|' 
G>(x) = x\/ x"^ — 1 — arccosh(rE) for x > 1, 



G<(x) = xy/l — x"^ — arccos(a;) for |x| < 1, 



where 



u = aq cos 6„, 

Vr. 



U={1- a'^){q^ - 1/2) 

2^ + bJ 



Xl 



X2 



Vc{q + au cos 6q) 

2fi — uj 
Vc{q + au cos 6q) 

_|_ 2^ + bJ 



VrU 



{q + OLV cos dq) lb 



a\ sm ( 

w 



1 



72 — v'^\ 



(2^ + 1 



U 



(B-2c) 

(B-3) 

(B-4) 
(B-5) 

(B-6) 
(B-7) 
(B-8) 
(B-9) 

(B-IO) 
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\q + av cos 9q\ it 



= a — cos I 



l-a2 
2 



a I sin ( 
W 



+ 



1 



4q;^ cos ( 



l-a2 



+ 



2a 



1 - a2 / ' 
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(B-12) 



a cos Pg + 



(+)1 



II /i *1 — a2 y \ fi J 1 — q2 



2 4a 
+ 



cos ( 



+ 



2a 



1 - a" 



(B-13) 



Appendix C: Expression of the Real Part 

The real part is calculated using eqs. (20) and (21) in which the semianalytical calculation 
is performed by dividing the part into four regions, i.e., (1)0 < Vcq/fJ^ < 1/(1 + a|cos0g|), 
(11)1/(1 + a|cos0g|) < Vcq/fJ. < 1/(1 - a|cos0q|), (111)1/(1 - a| cos 6*^1) < Vcq/fJ. < 



(C-l) 



2/(Vl - a2^1 -a2 cos2 9q), and (IV) 2/(Vl - "2^1 - a2 cos2 Oq) < Vcq/f^ 

Ren(g,0q,a;) = n^B+ ^ n'^, 

C=i,ii,ni,iv 

where the first term denotes the the contribution of Hgg, while IIj, ...,njy are the semianalyt 
ical expressions including the integral. 
The final result is given by 

rr' ff \nf 2 2m q{uj - ujBiOq)} 

^3B =f{(l,^)Q{^ - g )log 



'{u: - Wb(^<j)} + g2 _ 1.2 + ^(j,2 _ ^2) I (^g(0^) _ aqCOsOqf - g2| 

q{u: + ujb{t^ + Oq)] 



{uj + WB(vr + 9q)] + q^ - u'^ + J {v"^ - q^) { (wb(^ + Oq) + aqcos Oqf - q^^ 



+ f{q,i^)&{q —V )|^arcsin 



f . u{uj - ujB{Oq)] + q^ - y"^ . v{u: + u)b{t^ + 6q)] + q 
{ arcsm , — h arcsm , ; -; 



2 1.2 



g|cj - OJB{0q 



ni =|l+A^cs(0q),O] +I^^[up^{eq),UB{eq)\ + -lTaIl{eq,u)\0g 

+ /f^ ^cs(7r - Oq), O] + /2B [t^A(7r - 6*^), tJB(vr - Oq)] 

X e 



q\u + ujb{t^ + Oq) 
(C-2a) 

U: - UJresiOq) U} - UJ^iOq) 



UJ UJ - UJA{Oq) 



1 + a cos ( 



n'n 



I^j^[uJA{Oq),0] + I^j^[uJrcs{Oq),UJA{Oq)] + ^Bi^rcsiO g) , UJb {0 q)] + -lmU{Oq, u)log 
+ I-j^ ^cs(^ - Oq), 0] + [uJa{tT - Oq), UBiTT - Oq)] |© (| " ^g) 



(C-2b) 

UJ - UJB{Oq) 



UJ 
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+ I^js^[u}A{7r - 9q),0] + l2AKcs(vr - Oq),UJA{T^ " Oq)] + [^rcs{T^ " 6lq),WB(vr - 6q)] j© (Oq 



IT 

2 



X G 



1 — Q COS ( 



Vr.Q \ & \ V, 



1 + a cos ( 



=|AA[^A(^q),0] +4A['^rcs{0q),OJA{0q)] + Kcs (^q) , ^^B (^q)] + -Im 0(0^, w) log 

+ [a;A(7r - 6*^), O] + I2A [wres(vr - 6',), WA(vr - 61,)] + Ige Kes(vr - ^g), ^jbItt - 6*5)] | 
2/x 



(C-2c) 

W - UJBiOq 



UJ 



X e 



^(l-Q2)(l-a2 cos2 9g) 



Vcq]Q{ Vc 



1 — a\ cos ( 



(C-2d) 



n'lv =y^A['^resi9q),UJ+{9q)] + /2+B [^rcs (fq) , (^q)] + -Im 0(0^, w) log 
+ /2AKcs(7r - 9q),U}+{TT - 9q)] + ^cs(7r - 9q),U}B{TT - 9q)] | 



UJ - UJB{9q) 



UJ - UJ+{9q) 



2^1 



V^(l-Q2)(l-a2 C0S2 Oq) ) 

Here, we denote the functions defined by the integrals 



Aa/2a["'^] = f 
/2+B [a, b] 



b — aq cos Oq r\ \ / r\ \ 

arccos ; -f^iA/2A(9 COS y + aq COS 9q) — sin ylmll[9q,uj) 



dy 



a — aq cos oq 



q{cos y + a cos 9q) — oj 



arccosh 



b—aq cos Oq 



dy 



77- / , a — aqcos&n 

' ' J arccosh — 



H2B{qcoshy + aq COS 9g) — sinhylmn(6'q, 
g(cosh y + a cos 9q) — uj 



9 

h-\-aq cos Oq 

arccos - - HiA/2A[q COS y - aq COS 9q) 

dy 



^ a + aq COB Oq q{cOS y — a COS 9q) + UJ 



arccosh 



b+aq cos 6 



dy 



H2B (q cosh y — aq cos 9q) 



71" Jarccosh^i^i^^ g(cosh y - acos9q)+uj ' 



(C-2e) 

(C-3a) 
(C-3b) 
(C-3c) 
(C-3d) 



where 
HiaH 

H2b{^^) 



167r 



167r 

q 

167r 



G> {xiiuj)) + sgn {xi{uj) - 1) G> (x^ (w)) - G> - sgn {x2{uj) - 1) G> {x^ (uj)) 

(C-4a) 

G> (x+(a;)) + sgn {xi{uj) - 1) G> (x^(a;))] , (C-4b) 
vr + sgn (x2(tj) - 1)G< (x^ (a;)) -sgn (2:2(0;) + 1)G< (-2:2(0;)) . (C-4c) 
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Fig. 7. Normalized imaginary part, lmn(0q, a;)i;^/(7sM on ttie plane of Vcq/ fJ- and w/^ for 9q = Q (a), 
= 7r/2 (b), and 9q = tt (c), where a = 0.8. The straight line including the origin represents the 
resonance frequency Wros. The dark (bright) area represents the large (small) lmll{9q,uj), while 
Imn(6'q,c<;) = in the white area corresponding to the regions of 3A and IB in Fig. 2. The 
dash-dotted line (dotted line) denotes the boundary for the excitation process at which the cusp 
emerges (is absent). 
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Fig. 8. (Color online) Real part Kell{q,9q,ijj) as a function oiu/fi with a fixed 6q = tt/2, for qVc/fJ- 
= 2 (a), l/2(b), and 4(c), which correspond to Figs. 3 and 5(a) and 5(b). The global feature with 
arbitrary q is seen for uj > Wros- The divergence of the real part for w — > Wros — is obtained for 
q > 1/(1 + acos0,j). The dotted line represents the isotropic case of graphene. 
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Fig. 9. (Color online) Normalized interband optical conductivity Re(j(0g, cj)2/i/7re^ as a function of 
(jO / ^ for 6q = 0, 7r/4 and 7r/2. The contribution comes from only the interband excitation. Note 
that Jiea{9q) = Re(T(7r — 6q). The dotted line followed by a jump represents the isotropic case of 
graphene. 



(b) k 




Fig. 10. (Color online) (a) Particle-hole excitation with zero momentum transfer from the lower cone 
to the upper cone where three kinds of ellipses lead to the energy with u = 2/i/ (1 + a), 2/^ and 2 
/i,/(l — a), (b) The three kinds of ellipses of (a) are projected on the k^-ky plane. The excitation 
is allowed only in the region which is located outside the Fermi surface (solid ellipse). 
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11. (Color online) Plasma frequency as a function of VcQ/ H for 9q — 7r/2, 7r/4, 37r/4, 0, tt (from 
the top to the bottom). 
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12. (Color online) Real part at a; = as a function of Vcq/ /i with fixed 9q ~ 0, 7r/4, tt/2, where 
the lowest line denotes the result of the isotropic case of graphene. 
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